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We analyze the electromagnetic response of a system of charged bosons cou- 
pled to a Chern-Simons gauge field. Path integral techniques are used to obtain 
an effective action for the particle density of the system dressed with quantum 
fluctuations of the CS gauge field. From the action thus obtained we compute 
the U{\) current of the theory for an arbitrary electromagnetic external field. 
For the particular case of a homogeneous external magnetic field, we show that 
the quantization of the transverse conductivity is exact, even in the presence of 
an arbitrary impurity distribution. The relevance of edge states in this context 
is analyzed. The propagator of density fluctuations is computed, and an effective 
action for the matter density in the presence of a vortex excitation is suggested. 
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1 Introduction 



Quantum Field Theory (QFT) techniques applied to Condensed Matter systems have 
been extensively used in the recent past JTJ f2|. Specially, path integral methods have 
been applied successfully over the years in many problems of Condensed Matter physics 
such as superconductivity, superfluidity, plasmas, cristalization Q and quantum Hall 
effect (QHE) M. The essence of this approach is the evaluation of an effective partition 
function, obtained from the original one by integration over some degrees of freedom. 
In general, this technique leads to the problem of solving a fermion or boson functional 
determinant, to be evaluated by some approximation scheme. An application to the 
study of two-body interactions in the QHE can be found in ref. 0. Also, an interesting 
application to quasi one-dimensional interfaces was developed in ref. |J. 

In the last years, models involving matter fields coupled to a £7(1) Chern-Simons 
gauge field were analyzed, not only as QFT models 0, but also as fenomenological 
Condensed Matter models||. The most interesting feature of these models is that the 
Chern-Simons field attaches quantum fluxes to the particle positions (or to the particle 
density distribution) thus changing the statistics of the original fields. In particular, if 
we attach an even number of quantum fluxes to a matter field, its statistics is preserved 
(bosonic or fermionic). On the other hand, if we attach an odd number of fluxes, the 
statistic of the relevant degrees of freedom changes from bosonic to fermionic and 
vice- versa. This means that a system of fermions can be regarded as either a system 
of bosons with an odd number of attached fluxes or a system of (another kind of) 
fermions plus an even number of fluxes attached to them. The first possibility has been 
explored in the construction of the Chern-Simons-Landau-Ginzburg (CSLG) theory 
for the Fractional Quantum Hall Effect (FQHE)[[|; the second one has led to the 
"composite fermion" model for the FQHE [|I0) |T| . 

The electromagnetic response of the CSLG theory was calculated in ref. pi. In 
that paper, the Hall conductivity is evaluated in the saddle point approximation plus 
gaussian fluctuations. The essence of this approximation scheme can be described 
as follows: first, the path integral of charged bosons in the presence of an external 
gauge field is evaluated in the gaussian approximation, yielding a quadratic effective 
action for the external gauge field; then, the Chern-Simons coupling is turned on, and 
the integral over the statistical gauge field is computed. In this way, an interesting 
relation between bosonic superconductivity and incompressibility of fermionic degrees 
of freedom is established. 

Several authors argue that the quantization of the transverse conductivity should 
be exact due to the topological character of the Chern-Simons term. Nevertheless, to 
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the best of our knowledge, a complete explicit calculation, involving random impurities, 
is still lacking. In the case of the Integer Quantum Hall Effect (IQHE) (where electron- 
electron interaction is not considered), it was shown that the Hall conductivity is not 
affected by a delta-type impurity fT2l Hl3f . However, the study of a general type of 
impurities combined with a Chern-Simons gauge coupling is a very difficult task and 
no exact results are available. 

Another aspect of these models is the study of density fluctuations around the 
ground state. At the classical level, topological [[Uj][]n| [|16| and non-topological |L7 



solutions to the equations of motion were found. However, it would be interesting 
to study if the quantum fluctuations of the gauge field do really modify the density 
profiles. 

Motivated by these ideas, we will consider in this paper the electromagnetic response 
of non-relativistic bosons coupled to a Chern-Simons Gauge Field. To this aim, the 
path integral technique is employed, so that all degrees of freedom of the model but 
the particle density are integrated out; hence, the resulting partition function — a 
functional integral over the charge density alone — takes into account the quantum 
fluctuations of the gauge field. (The charge density is, perhaps, the simplest of gauge 
invariant objects.) In particular, the Chern-Simons constraint (that attaches fluxes to 
particles) is implemented exactly, rather than being artificially imposed by means of 
some approximation scheme. 

This technique allows us to clarify some aspects of this system. First, we can write 
an explicit expression for the U(l) current of the model submitted to an arbitrary 
electromagnetic field. Specializing this general expression for the current to the case 
of a static and homogeneous external magnetic field, we can explicitly show that the 
quantization of the transverse conductivity is exact, even in the presence of any type of 



impurities. Thus, we generalize Prange's result [|12j||T3[ to the FQHE and to any kind 
of impurities in the context of the CSLG theory. 

This paper is organized as follows: In we deduce the effective action for the 
bosons coupled to a Chern-Simons field plus an arbitrary electromagnetic field, in the 
density representation. In §|3| we evaluate the U(l) current of the model. We show the 
exact quantization of the Hall conductivity and the role played by the edge states. In 
§|] we analyze the dynamics of density fluctuations around the ground state. Finally 
we discuss our results and present our conclusions in §[5|. 
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2 The Effective Action for Bosons Coupled to a 
Chern-Simons Gauge Field in the Density Rep- 
resentation. 



Let us begin by considering the Euclidean action for non-relativistic bosons coupled to 
a U(l) Chern- Simons gauge field a^x) and an arbitrary external electromagnetic field 

M x )> 



S = f d 2 xdr I 4>* (d T + i[Ao + a ] -ifi)<j>-— (- V - e(A + a) 
J [ m \i 

- X -j d 2 xd 2 ydr (0*(z)0(z) - p) V(x - y) (<F(y)<Ky) - p) + 

- \(^)J d 2 xdr e^a,d v a p 



+ 



(1) 



where V(x — y) is an arbitrary two-body potential, and the last term is the well-known 
Chern-Simons term. Due to this term, the classical equation of motion related to ag 



7T. 



p(x) = 4>*4> = ^-V x a 



(2) 



is a pure constraint that attaches | fluxes to the particle density. In order to represent 
fermionic degrees of freedom from the bare bosonic fields, we must choose 9 = (2n+l)7r, 
with integer n. 

The partition function is 



Z(A ,A l )= fv<irV4>Va lt G F {4>,a ll ) e - s ^*' a ^ 



(3) 



where Gf{4>, ct M ) is a gauge fixing functional, necessary to avoid overcounting of physical 
states in the functional integral. 

The presence of the Chern-Simons field allows us to extract the external electro- 
magnetic field A^ from the kinetic operator of the action. Making the following change 
of variables: 



Va, 



<V - 
Va,. 



A,. 



(4) 
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we obtain 



S 



d 2 xdr < <p* (d T + ia — in) <\> 

m 



-V + a 

% 



i j d 2 xd 2 ydr (<P*(x)<l>{x) - p) V{x - y) (<f>*{y)<f>(y) - p) 



l 7T 

29 

7T 



d 2 xdr e i0j 'a,-a 



7T 



d 2 xdr oo(V x a + B) 



— 2 







(5) 



where the matter field couples to the dynamical gauge field only. Here we use E and B 
to denote the external electric and magnetic field, respectively. (In what follows, greek 
indices run from to 2 and latin indices run from 1 to 2.) 

Going over the density representation, it is useful to perform another change of 
variables, 



with trivial Jacobian 



Replacing in (|) we have 



<j)(x) = pi/*( x ) e ie{x) 



(6) 
(7) 



S 



d xdr < -d T p + ip(d T 9 + a )> - - — < — Vp ■ Vp + p 

2m I 4p 



1 f 1 



V0- a 



d xd ydr (p(x) — p) V(x — y) (p(y) — p) + p / d xdrp(x) 



1 

2 

-— / d 2 xdr e ■'aidoaj — i— / d 2 xdr ao(V xa + B) 
2 9 J 9 J 
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The idea is to integrate out the gauge fields ao(x), aj(x) and the phase 9(x), in order 
to obtain an effective action in terms of the density p(x) only. At this point, it is useful 
to discuss the gauge fixing functional Gf{4>, ffy)- 

There is great freedom in choosing Gp{4>i a n), and the partition function Z(A , AA 
must be gauge invariant, i.e., it must not depend on any particular choice. Notice that 
in the "polar coordinates" we are using, a gauge transformation reads 
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6(x) — > 9(x)+A(x) (9) 
a^x) — > a^x) + d^A(x) . (10) 

So, we can always choose A(x) = — 6{x) to have 9{x) = in the action. This corre- 
sponds to choosing the fixing functional to be 

G F (9,a^)=S(9(x)) . (11) 

We must prove that this is an accessible gauge, and that it does completely fix the 
gauge. It is not difficult to show H that this gauge is equivalent to Coulomb's gauge 



G F (6,a^ = 8[V-a). (12) 

Demonstrating such equivalence consists of showing the existence of a continuous gauge 
transformation that leads from condition ( TTf) to ([12]). Decomposing the field Oj into a 
transversal af plus a longitudinal component af, it is possible to rewrite the integration 
measure in the Coulomb's gauge in the following form, 

J Va{DpD0 8 (V • o) . . . = / Va^VpVQ ... (13) 

as the Coulomb's gauge sets the longitudinal component of the gauge field to zero. Let 
us now perform a gauge transformation 

di(x) — > a,i(x) - di0(x) . (14) 

This transformation removes the phase of 4>(x) from the action but adds a longitudinal 
component af (x) = —di6{x) to af{x). Then, we can replace T>9 with T>a L {x) in ([l3| ) 
(with a trivial Jacobian) obtaining 

(13) = J Va L Va L (x)Vp... = J Va{DpV6 5 (6(x)) . . . (15) 



In this way, we show that ( |Tl~D and (|T2| ) are equivalent. 

Now, we are ready to begin with the integration. We see that equation @ has 
only linear terms in a (x). Integration is trivially carried out, yielding the well-known 
Chern-Simons constraint 

/ Va e~ l f d2xdT «o(p-| [Vxff+B]) = 5(jJ _ X a + B]) . (16) 
J 9 
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The quantization of this model in the framework of constrained systems has been 
considered by several authors; in particular, we find the gauge invariant approach of 
reference JHJ] very instructive. However, in the present paper, we prefer to follow a 
functional approach, i.e., we will integrate over all configurations of p and a, that satisfy 

p( x ) = ^[Vxa + B). (17) 

Note that the attachment of magnetic fluxes to the particle density p(x) is not a 
classical effect. It remains in all the integration over p and a^. A very interesting 
physical realization of this phenomenon can be found in [IS|. Thus, using (|T^) and 
(|TTD we find for the partition function 

Z{Ao,Ai) = JvpVa5(p-^[V xa + B]\ (18) 



with 



Sf = — - — J d 2 xdr I — Vp ■ Vp + p\a\ 2 1 + p J d 2 xdrp(x) 



4p 

d 2 xd 2 ydr (p(x) - p) V(x - y) (p(y) - p) 

-— I d 2 xdr e^didodj —i-zl d 2 xdr e^aiEj 
2 9 J J 

J d 2 xdr e^ p A^d u A p (19) 
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Clearly, Sf is not gauge invariant, as we have fixed the gauge to perform the integrals. 
Of course, after the integration, this symmetry must be restored. Another observation 
about Sf is that the density p(x) has no independent dynamics; rather, all of its 
dynamics depends on its coupling to the gauge field. This is an important fact in the 
study of the dynamic of density fluctuations. We will return to this point in section 

m ' 

The next step of our development is to integrate out the field aj, to obtain an effec- 
tive action only for p(x). In order to do that, it is useful to decompose the bidimensional 
vector ai into a longitudinal plus a transversal part 

a, = ditp + eijdjT] . (20) 
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This is a linear transformation. So, 

Va{Da 2 — ► VrfDip (21) 

This decomposition is equivalent to describing a vector field through its rotational and 
divergence, since 

V-a = VV (22) 
V x a = -V 2 r] (23) 

Replacing fl2"UD in fllpp we find 

= / d 2 xdr I — Vp • Vp \ + u f d 2 xdrp(x 

2m J 4p J 



I I , , . If J_ 

4p 

d 2 xd 2 ydr (p(x) - p) V(a; - y) (p(y) - p) 



- — / d 2 x<ir p(x) fv<p ■ V<p + V77 • V77) - ( d 2 xdr ipd V 2 i] 
+ i*- J d 2 xdr <p(V xE)+ r/(V ■ E) 



~J d 2 xdr eTPA&Ap (24) 

This action is quadratic in the field <p. Moreover, the S — functional in ( [TED depends 
only on 77. Therefore, we can integrate over <p using the relation 

J)<p e -J dx \f{x)dip{x)+.j{x)ip{x) _ Def-yi^o^i J dxd v J{ x )o~ x {x-y)J{y) ^ ^5) 
thus obtaining 



Z(A ,A l ) = JvpV V 5(p-^[B-V 2 V ]) e- 3 '? (26) 



with 



Sp = — — / d 2 xdr I — Vp ■ Vpl + p / d 2 xdrp(x) - -Tr ln(V • —V) 
2m J I4p J j 2 m 



d xd ydr (p(x) - p) V(x - y) (p(y) - p) 
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/ d 2 xdr p(x)V?7 • V?7 

2m J 



in ix 2 /* — > — * 1 — » — » 

+ T^2 / (d V 2 V (x) - V x E(x)) — (QoV^^) - V x E{y)) 

4 0* J V ■ (pV) 

+ %- Q J d 2 xdr n(y ■ E) - i~ y d 2 xrfr g^^AA, (27) 

where we have used the fact that InDetO = Tr In O. 

Finally, we must integrate over the transversal component of dj, say ^(x). This 
is not a difficult task because r)(x) is in the argument of a 5 — functional (see eq. 
(|I6D). From flTTD we see that the support of the 5 — functional consists in all the field 
configurations that satisfy 

V 2 r? = B - -p (28) 

7T 

We can invert this linear differential equation by using 

v{x) = J_f d 2 y G{x _ y)(p(y) _ * B{y)) (29) 
where G(x — y) is the Green function for the Laplacian operator, 

V 2 G(x -y) = 5(x - y) (30) 
The well-known property of the 8 — function 

*(/W) = E^Jp. /(^) = 0. (31) 
can be extended to the 5 — functional, yielding 

S (p - ^(B - V 2 ^)) = ^5 (v(x) + 9 ~Jd 2 y G(x - y) [p{y) - j;B{y)fj (32) 

The factor |V 2 5| is field- independent and can be absorbed in the global normalization 
constant of the partition function. Using (|3~2|) , it is now simple to functionally integrate 
out the field r], obtaining 

Z(A , Ai) = e^I d2xdT ^"^a p J Vp e -s eff (P,S,B) (33) 



9 



with 



- ^ y d 2 xd 2 ydr (p(x) - p) V(a; - y) (p(y) - p) - ^Tr ln(V ■ — x 
1 r d 2 xd 2 ydr [p{ x )- 7 LB{x))F{x-y)[p{y)-^B{y)) 



2 



\ 9 f -> -> 
J d 2 xd 2 yd 2 zdr VG(x - y) ■ VG(x - 2) x 

x (p(x) - |s(ar)) (p(</) - ^fi(y)) (p(z) - 

y d 2 xd 2 ydr (p(x) - jB(x)\ G(x - y)V ■ E(y) (34) 



9 1 r 

F(x-y) = -—J d 2 z B{z)VG{z - x) ■ VG{z - y) (35) 



where 



Equations (|33D and ( p4|) are the main results of this section. S'e// is the effective 
action for the Chern-Simons Bosons in the density representation, interacting with an 
arbitrary electromagnetic field. This is a gauge invariant action, as it only depends 
on p(x), E and B. The coupling of the charge density to the electromagnetic field 
is very peculiar. The magnetic field acts as a "background density", entering the 
action only through terms of the form (p{x) — |S(x)V The electric field couples to 

the density only through its divergence V • E, which is obviously proportional to the 
external charges (impurities, for example). So the last term of (0) is simply the two- 
dimensional Coulomb energy between the external charges and the charge fluctuations 
over the magnetic background. In this term we absorbed the chemical potential as it 
can be simulated by a uniform background charge density. Note that, in particular, a 
homogeneous electric field do not couple to the matter field. 

The coupling of the gauge field with the bosonic matter field had two main con- 
sequences in the particle density of the system. The integration over the longitudinal 
gauge field produced a non-local dynamical term for p(x) (see first term of (|34D ) and a 
non-local interaction term given by Tr ln(V • — pV). The integration over the transver- 



sal gauge field induced a two-body interaction (third term of (|34|)) and a three-body 
interaction also. The induced two-body interaction is the basic characteristic of the 
fluctuation of the transversal degrees of freedom of Chern-Simons field and, for the 
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particular case of B(z) = B = const., 



F(x-y) = ---G(x-y) , (36) 
7r m 

we obtain the well-known logarithmic two-body interaction, extensively explored in the 
past. This interaction is responsible for opening a gap in the spectrum of excitations 
over the ground state of the system. Moreover, the longitudinal fluctuations are non- 
trivial. In particular, the non-local dynamical term will modify the propagator of 
the density fluctuations. Also, we note that a Chern-Simons term for the external 
electromagnetic field has been factored out from the action. In the next section we 
will take advantage of these properties to analyze the electromagnetic response of this 
system. 



3 Electromagnetic Response 



This section is devoted to the calculation of the electromagnetic response of the theory 
described in the last section. 

Due to gauge symmetry, the conserved current is given by 

<J ^ >= 6XM HZ) ■ (37) 

Explicitly evaluating the functional derivatives, we find using ([33|) and fl34l): 

<Jo(x)> = i lB + a i < SS '"}l?> B) )> (38) 

<Hx)> = -£e ijEj + d < SS '»& B) ) > < SS eff (E, B) ) > 

U orlji on 

— * 

As S e ff is gauge invariant, it depends on A{ only through E and B. This implies the 
existence of a topological current, automatically conserved (last term of ([39|)). We will 
see that this topological current is responsible for the exact quantization of the Hall 
conductivity, even in the presence of impurities, a fact closely related to the gauge 
invariance and edge state excitations. 

We can easily evaluate expression ( [38|) , obtaining the obvious result 

< i J >=< p(x, t) > (40) 
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The calculation of is direct, but the result is less trivial: 

<J l > = -z^ey - dj J d 2 y G(x - y)V ■ E^j 

+ i J d 2 y diG(x - y)d < p(y) - ^B(y) > 

+ J?{x) (41) 

where 

jj{x) = —^Aix) (42) 

and 

A(x) = J d 2 zd 2 y VG(x - y) ■ V y G(y - z) < (p(y) - ^B{y)){p{z) - ^B(z)) > (43) 

This is the first result of this section. The current (|4T|) is the response of the 
system to an arbitrary electromagnetic field. Before we go over to the calculation, it is 
necessary at this point to check the conservation of the current ([4~ID . So, let us evaluate 
the divergence of < Jj >. The last and the second terms of the first line of (|4"ID are 
automatically conserved because both of them are topological terms. That means, they 
are of the form €ijdjf(x), with arbitrary /. The other part of the current needs some 
attention. From ( [41]) we have 

d t < Ji >= -i^eijdiEj +tjd 2 y V 2 G(x - y)d < p(y) - ^B(y) > (44) 

Using ( |30D we can integrate the last term obtaining 

di < Ji >= id < p{y) > -i~ (VxE + d B) (45) 

From fllof) (remembering that the external electromagnetic field satisfies Faraday's 
Law), we obtain the current conservation law (in Euclidean space) 

di < Ji> +d < J >= (46) 

For an arbitrary electromagnetic field, expression ( fill) can turn out to be quite a 
complicated one. However, it is useful for the study of the structure of the current 
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for specific configurations of the electromagnetic field. For example, let us analyze the 
case of a static and uniform magnetic field, and a static but arbitrary electric field. 
This case is the relevant one in the study of the Quantum Hall Effect. In this case, 
equation ([41]) reduces to 

< Ji > = (Ej - dj J d 2 y G{x - y)V ■ E^j 

+ J?(x) (47) 

The first line of (47) does not depend on any specific detail of the system. Clearly, 
the dynamical properties of the system (two-body interactions, etc.) are all contained 
in the topological current Jj(x). So, let us analyze this current more carefully. From 
we have 

A(x) = J d 2 zd 2 y VG(x - y) ■ V y G(y - z) < 5p(y)5p(z) > (48) 



where the expectation value is evaluated with the partition function (|33"D, and 5p is the 
density fluctuation around the magnetic background. The action S e ff only depends 
on E through its divergence. So, in a clean sample (without external charges), A (re) 
does not depend on the electric field. Moreover, we have that < p > is constant due 
to translation invariance and < p >= p fixes charge neutrality. Also, we saw that the 
effect of the quantum fluctuations of the Chern-Simons field was that of inducing a 
new two-body and three-body interactions with neutralizing background ^B. In order 
to define the density fluctuations Sp of equation (ESI) correctly, we must have 



p = - Q B (49) 

thus implying a Landau filling factor v = |. In this case the two-point correlation 
function of density fluctuations < 5p(x)5p(y) > depends on \x — y\. In the case of 
a sample with impurities, we must calculate the average of the correlation function 
with some weight factor (gaussian, for example). Thus, upon averaging the random 
impurities we restore the translation invariance and 



< 5p(x)5p(y) > = a(x- y) (50) 
(the bar means "average over random impurities"). In any case, we can write for A, 

A(x) = J d 2 zd 2 y VG(x - y) ■ V y G(y - z)a{y - z) 

' d 2 y VG(x -y)-V y J d 2 zG{y - z)a{y - z) 

(51) 
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In an infinite plane, the last integral of ( |5"ID is constant and therefore A(x) = 0, for 
v = tt/6. So, in an infinite system the topological current jj given by (|42| ) is zero. 

In a finite plane, jj is no longer zero, but we will show that it does not contribute 
to the total current. To see this, let us write the total current due to the topological 
density current as 

I T = [ jf dS=—-[ e^Aix) dS (52) 
Jd mv: Jd 

The domain of integration D is simply the surface of the sample. Supose that D is 
limited by a border dD. In this case we can represent the function A as A(x)Q(D), 
where 0(D) is one in D and zero otherwise. So, 

dj(A(x)e(D)) = d j A(x)Q(D) + A(x)djQ{D) 

= d j A(x)Q(D) -n j A(x)5(dD) (53) 

where 5(dD) is the Dirac delta-functional with support in dD, and hj is the j-th 
component of the versor contained in the plane, normal to dD and external to D. 
Hence, 

i e 



if = — - || e ij d j A(x)Q(D) dS - J e ij n j A(x)8(dD)\ 

dltiA(x) - J e ij n j A(x)5(dD) ds\ (54) 



i e 

where U = e tangent versor to dD. But, 

/ e ij n j A(x)5(dD) ds = <f dli^x) (55) 

J J 3D 

which implies that If — 0. 

The r.h.s of ([55]) represents an edge current, responsible for the cancellation of the 
topological total current. It means that, due to the edge currents, the main contribution 
to the total current comes only from the first line of fl4"T|). The relevance of edge states in 



the QHE was first stressed in |T9[ and then, its theory was developed in refs. |20[ ||2T| |22] . 

It is interesting to note that, for v ^ there are two different "uniform back- 
grounds", say, p and (it/6)B. In this case, it is not possible for the mean value of the 
density being constant. The translation invariance is broken and the density fluctua- 
tions no longer depend simply on \x — y\, but show a non-trivial dependence on x and 



y. This shows that for v ^ tc/8, the current will not be given by the first line of ([4lD 
but the topological density current Jf will contribute in a non-trivial way. 
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Let us analyze now the main part of (pLTD, namely 

< Ji >= -i|ey ^ - a,- 1 d 2 y G(z - y)V ■ E^j (56) 

Its second term has support only in the region p(x) 7^ and depends on the divergence 
of E, only. Thus, in a clean sample, V • E — and the current is the well-known Hall 
current (turning back to usual units) 

< Ji >= ^~^E S (57) 

In a real sample, where impurities are present, V-£ 7^ 0. In this case we can decompose 
the electric field E in two parts 

E = E' + E imp (58) 

where V • E' = and V • E imp = p e , with p e the external charge due to impurities. 
With these definitions we can rewrite flSTf ) in the following form: 

< J 4 >= -i^eyE's - ijey [E^ - d s J d 2 y G(x - y)V ■ E imp ) (59) 

It is simple to show that 

Ef v - dj J d 2 y G(x - y) V • E imp = 0. (60) 

The second term of the l.h.s of the previous equation can be interpreted as the electric 
bidimensional field created by a density charge p e = V • E imp . 
Thus, from (|59"D we have (in usual units) 



< J - >= stt^ (61) 



So, the quantization of the Hall conductivity for v = 71/6 is exact in this model and 
does not depend on any impurity distribution. In other words, the transverse current 
only "sees" divergenceless fields. System dynamics and random impurities only affect 
the edge states of the system, and this the reason for the exact quantization of the 
Hall conductance. This result generalize Prange's result for a 5 — type impurity, 
to the case of a general impurity distribution in the context of the CSLG theory for 
the FQHE. 
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4 Dynamics of Density Fluctuations 



In the last section we calculated the electromagnetic current and, from the resulting 
expression, we were able to understand the role of the Chern-Simons term, of the edge 
states — that are responsible for the exact quantization of the conductivity — and the 
role of impurities, too. However, for a complete understanding of the electromagnetic 
response, it is necessary to evaluate the polarization tensor 

, . 5 2 lnZ(A) , , 

n ^ = JMyM) < 62) 

The calculation of the functional derivatives is straightforward and the result is formally 
expressed in terms of mean values of the density fluctuation < 5p(x)5p(y) . . . >. In 
particular, the expression for II 00 is very simple 

U 00 (x,y)=<6p(x)5p(y)> (63) 

were Sp(x) = p(x)— < p(x) > 

In order to explicitly calculate this type of objects we need to evaluate the effective 
action for the density fluctuation Sp. The action Q3~4] ) is the appropriate one for the 
study of the dynamics of density fluctuations around the ground state of the model. 
The structure of the ground state depends on the specific configuration of the external 
electromagnetic field. We will specialize our present analysis to the case of a static and 
homogeneous field, the relevant case to the study of the QHE. The ground state for 
v = tt/6 is homogeneous, and is given by 

< p(x) >=P = ^B (64) 



so that p(x) = p + 5p(x). The aim of the present section is to obtain from (|34|) an 
effective action for 8p{x) and to extract from it the Feynman's rules to calculate mean 
values. We need to perform a Taylor expansion of the action around the configuration 
p(x) = p. This expansion is straightforward, except for two terms originated from 
integrating over the longitudinal component of the gauge field. These terms are the 

kinetic term d^p j ^ j d^p, and the determinant obtained when integrating over tp, 

namely, Tr ln(V • — V). Let us analyze these two terms carefully. We can write the 
kinetic term of fl34]) in the following way: 



IK 



d 3 zid 3 z 2 Sp(zi)}C(z 1 ,Z2)Sp(z2) (65) 
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(66) 



where 

V 21 • (p{ Zl )V Zl K{z x , z 2 )) = 8{z 1 -z 2 ) 

(5p means the derivative of 5p with respect to time). We have also used that the ground 
state density is static. Note that K depends on p through the implicit equation 
Thus, we can expand the kernel K, in powers of 5p in the following form, 



From 



^Oi, z 2 ) = K.(z u z 2 ) 
we see that 



p=p 



+ d 6 x 



, 3 5K(z 1 ,z 2 ) 



5p{x) 



Sp(x) + 



p=p 



K.(zt, z 2 ) 



\p=p 



-G(zi, z 2 ) 
P 



(67) 



(6J 



where G(zi,z 2 ) is the Green's function of the Laplacian (see eq. (|30|) ). Also, we can 
show (see Appendix) that 



5K{zx,z 2 ) 



5p(x) 



p=p 



; VG( Zl ,x) ■ VG(x,z 2 ) 



(69) 



Introducing (|68|) and (^9]) into (|65|) we find 
1 



IK 



+ 



+ 



P 



P 2 



d 3 zid 3 z 2 5p(zi)G(zi ) z 2 )8p(z 2 ) 
d 3 zid 3 z 2 d 3 x VG(zi,x) ■ VG(x, z 2 )5p(z 1 )8p(z 2 )5p(x) 



The next interesting term to be considered is 



which can be rewritten as 



J de ; = Trln(V-^V), 

m 



U = Trln(^V 2 + V^V) 

m m 



(70) 



(71) 



(72) 



Observe that this term is due to the interaction of the density field with the longitudinal 
component of the gauge field. This process is governed by low-momentum transfer, so 
the logarithmic term with two derivatives can be neglected, and we obtain 



r de ^Trln(^V 2 )+Trln(l + ^) 
m p 



(73) 
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The first term in this expression is a constant and can be absorbed into the global 
renormalization factor of the partition function. The other term can be expanded in 
powers of Sp, yielding 

The linear term can be absorbed by a suitable renormalization of the chemical 
potential; the only effect of the quadratic term is redefining the two-body potential 
V(x,y) in the low-momentum limit. In the case of local potentials, this term may be 
used to redefine the coupling constant, but it is irrelevant in the case of long-distance 
interactions. 

Collecting all the terms, we obtain the effective action for the density fluctuations, 

TTl 1 C 

S = --Jd 3 xd 3 y5p(x)G(x-y)5p(y) 

+ \ J d 3 xd 3 y 5p(x) |^-V 2 5(x - y) - (v(x - y) - £(^) 2 G(x - y)) } 5p{y) 

+ j d 3 x V5p(x) ■ V5p(x)5p{x) - j d 3 x 5p 3 + ... 

1 2 f -> -> 

I d 3 xd 3 yd 3 z VG(x - y) ■ VG(x - z) 5p(x)5p(y)5p(z) 



+ 



2m 7r 2 
m 1 



TT? f — * — ♦ 

— — / d 3 xd 3 yd 3 z VG{x - z) ■ VG{z - y) 5p{x)5p{y)5p{z) 
2 p J 

d 2 xd 2 y 5p(x)G(x — y)V • E (75) 



As the fluctuation 5p is small, the quadratic term of the effective action may be 
regarded as defining the dynamics of the fluctuations, whereas the cubic terms stand- 
ing for perturbations to such dynamics. In this sense, the propagator of the density 
fluctuation in Fourier space may be written as 

< 5p(cu, k)Sp(-u;, -k) >= r — ^ — T (76) 

- 2 -{GS) +^(-k 2 ) + ^(f) 2 } 

This propagator has very interesting information encoded. For example, the compress- 
ibility of the system is defined as: 

k = lim n 00 (0, k) . (77) 
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From (|63| ) and the propagator ( [7TD we have for all the potentials with the general form 
V(-k 2 ) oc k a , 

« = lim ^ ~ = ( 78 ) 

showing that the system is in an incompressible state. 

The relation of this incompressibility to boson superconductivity was first shown 
in ref. 0. In the context of the present work, incompressibility depends on two facts. 
The first one is the existence of a gap in the spectrum of excitations. However, the gap 
existence is a necessary but not sufficient condition to determine the incompressibility 
of the ground state. The second fact is the presence of a k 2 factor in the numerator 
of the propagator. The gap in the spectrum is opened by the transverse gauge field 
fluctuations and its value 5 A = ^2-(^) 2 coincides with those found in refs. 0, |fLl 



and This gapfull excitation corresponds to the cyclotron mode excitation, and 

was identified as an inter-Landau level excitation. The k 2 factor in the numerator is 
induced by the longitudinal gauge field fluctuations that affect the dynamics of the 
density fluctuations. Remember that these fluctuations are completely equivalent to 
the phase fluctuations of the matter field. 

It is useful to write the propagator ([F(]) in the following form: 

< 5p(cu, k)6p(-u, -k) >= ^k 2 (— i— ) + 0(k 3 ) (79) 

V \U — UJk OJ + UkJ 

with 



f^V + A k 2 n-k 2 ) + 4(-)- (80) 



\ \2m 

This shows that our approach is consistent with Kohon's theorem ||23|| . This theorem 
asserts that in a planar translation-invariant system submitted to a perpendicular 
magnetic field. The density- density correlation function up to order k 2 must have a 
gap equal to cj c = — . This is an exact result, independent of the microscopic details of 
the system. Thus, a consistent calculus of radiative corrections must take into account 
that constraint at each order of the perturbation expansion, in addition to the vertices 



read from (^). Note that this result is valid only for two-body potentials that vanish 
as r — > oo. For example, a logaritmic potential V(— k 2 ) oc 1/k 2 modifies the gap, as 
can be confirmed from (|80|). 

The dispersion relation ([8(]) coincides with that calculated in |§ and in the 
low-momentum limit. For the Coulomb potential V(— k 2 ) oc l/|k| the dispersion is 
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linear, as it was first suggested in ref . |[24|| . Moreover, equation ( |50"D coincides with the 
dispersion relation of quasiparticles in an anyon superfluid [p^| . It is interesting to note 



that, in ref. [H] the model is written in terms of gauge invariant operators 

$(x) = Mx)^^ $(xY = e- l ^ x) Mx) 



where (p{x) is the longitudinal component of the gauge field given by (|20|). The theory 
is canonically quantized and then, the excitations of quasiparticles are calculated in the 
Bogoliubov approximation. In our functional approach, the gauge invariant degree of 
freedom is chosen to be the particle density p(x) and integration over all the remaining 
fields is carried out. Studying density fluctuations, we find that the dispersion relation 
of the cyclotron excitation is exactly the same as the former. This should mean that 
the term (k 2 /2m) 2 of ( |S0D is closely related to fluctuations of gauge invariant objects. 

It is well known that this model presents another type of excitations called topo- 
logical vortices. The typical asymptotic behavior of a vortex is (for \x\ — > oo) 

(j){x) = pi/2 e *r(*) (82) 



x 



cti(x) = da(x) = eij— (83) 



\x 



where 7 is the azimutal angle. Finite energy requirements conduce to flux quantization 

a-ds = 27v (84) 



and due to the Chern-Simons constraint the charge of the vortex is also quantized. 
Explicit solutions to the classical equations of motion were found in several references 
Hl5l | [jl6|| . Here, we would like to find an effective action for the particle density con- 



taining the vortex dressed with quantum fluctuations of the gauge field. The effective 
action (|75| ) has lost its topological information and contains no vortex. So far, we 
worked out all the functional integrations by implicitly fixing trivial boundary condi- 
tions on the fields. In order to obtain a density representation for the effective action 
in the presence of a vortex, though, we must fix non-trivial boundary conditions when 
performing the functional integrals. 

In ( |i~T|) we fixed the gauge by imposing S(6(x)), and we have shown that this is 
equivalent to choosing Coulomb's gauge. However, this is not the only possible choice. 
Supose we allow a singular gauge transformation, 

9{x) — > 6{x) + A(x) +7(x) (85) 
a^(x) — ► a^(x) + <9 M A(x) + d^{x) (86) 
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(7 is the azimutal angle) and choose A = —9. This transformation fixes the phase to the 
value 6 = 7 in the action, and induces a longitudinal gauge component af = d^B + c^y. 
Thus, we can change the integration variables from 6 to af. Retracing the steps that 



lead from ([TT|) to (0), we conclude that we can also choose the gauge fixing functional 
to be 

G F (8,a l )=5(6- 1 ) (87) 

This gauge fixing automatically establishes topological boundary conditions for the 
phase or, equivalently, for af . A similar technique was also used in P5 |, in the context 
of the abelian Higgs model. 

With these boundary conditions, the action in terms of density fluctuations takes 
the form 

S v (p) = S eff (p) + S^p) (88) 
where S e ff is the effective action with trivial boundary conditions given by (|34j) and 

S» = i / d 3 x {-\p\ V 7 | 2 + pa ■ V 7 ) ■ (89) 
Integrating over the gauge field, we obtain 

S 7 {p) = [ d 3 x p\Vj\ 2 + — Jey / d 3 xd 3 y d iP {x)G{x - y)j(x)d jP (y) (90) 

The first term of this equation comes from the integration of the longitudinal fluctuation 
of the gauge field, and can be interpreted as the interaction with an external charge 
density, as it can be rewritten in the form 

/ d 3 x p|V 7 | 2 = - / d 3 xd 3 y p(x)G(x - y)p e (y) (91) 
Am J m J 

where p e (x) = V 2 |V7| 2 = l/|x| 4 . 

The second term of (|90| ) comes from the integration of the transverse gauge field 
fluctuation; such fluctuation induces a new two-body interaction relating density gra- 
dients in orthogonal directions. Note that S" 7 is not translationally invariant because 
7(x) is not well defined at the origin (the nucleus of the vortex). Thus equation fl8~8|) , 
together with fl34f) and (0), is the effective action for the density charge of bosons 
coupled to Chern- Simons gauge field in the presence of a vortex centered at the origin 
of coordinates. It would be possible, upon minimizing this action, to obtain a vortex 
profile dressed with the quantum fluctuation of the gauge fields. 
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5 Discussions and Conclusions 



We have considered in this paper a system of non-relativistic bosons coupled to a 
Chern-Simons gauge field and an arbitrary external electromagnetic field. We have 
used path integral techniques to deduce an effective action in terms of the matter field 
density only. We found the coupling of the charge density with the electromagnetic 
field very interesting. The magnetic field acts as a "background density", as it only 
enters the action through terms of the form (p(x) — ^B(x)^j. The electric field only 

couples to the density through its divergence V • E, which in its turn is proportional 
to the external charges (impurities, for example). So, in the absence of impurities, the 
external electric field does not couple with the matter field. This peculiar fact is due to 
the Chern-Simons structure of the dynamical gauge field and it is the main reason for 
the exact quantization of the transverse conductivity. 

The coupling of the gauge field with bosonic matter had two main consequences 
on the density of the system. The integration over the longitudinal gauge field pro- 
duced a non-local dynamical term for p(x) and also a non-local interaction given by 
Tr ln(V • ^pV). The integration over the transversal gauge field induced a two and 
a three-body interaction terms. The induced two-body interaction is the basic char- 
acteristic of a fluctuation of the transversal degrees of freedom of the Chern-Simons 
field, and it is responsible for opening a gap in the spectrum of excitations. Moreover, 
the longitudinal fluctuations are also non-trivial. In particular the non-local dynam- 
ical term modifies the propagator of density fluctuations, leading, in addition to the 
appearence of a gap, to an incompressible ground state. Another interesting aspect of 
this action is that, apart from the interaction terms, a pure Chern-Simons term in the 
external electromagnetic field is factored out. This fact is important in the quantization 
of the transversal conductivity. 

Using these properties of the action, we were able to build an exact expression 
for the conserved U(l) current of the model. In particular, we have explicitly shown 
that, in the case of a homogeneous and static magnetic field, the quantization of the 
transverse conductivity is exact, even in the presence of any type of impurities. Thus, 
we generalized Prange's result [|HJ for the FQHE to any kind of impurity distribution, 
in the context of the CSLG theory. This development is based on the fact that all 
the microscopic dynamics of the model and the coupling with impurities are entirely 
contained in an automatically conserved current (topological current) that does not 
contribute to the total current. In the case of an infinite system, this topological density 
current is zero. In a finite system the edge states exactly cancel out the contribution of 
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this density to the total current. Moreover, we have shown that the transverse current 
only "sees" divergenceless field, thus canceling any contribution from external charges. 

We have computed the propagator of density fluctuations, using our effective action. 
The propagator shows a gap in the excitation espectrum that coincides with the inter- 
Landau level excitations calculated previously using semiclassiclal approximations ||. 



The dispersion relation of this excitation coincides with that calculated in ref. |14| by a 
gauge-invariant canonical formalism in the framework of anyon superconductivity. This 
observation suggests that the form of the dispersion relation is related to fluctuations 
of gauge invariant objects. Of course, in the low-momentum (long-distance) limit, the 
dispersion relation is linear (as noted in p^j ) and coincides with all the other approaches 
in this limit. 

We have also analyzed the structure of the effective action in the presence of a 
vortex excitation. The main difference in this case is that the integration over the 
longitudinal gauge field produces an interaction with an induced external density charge 
of the form p e = l/|x| 4 . The integration over the transversal gauge field induces a new 
two-body interaction between the gradients of the density in orthogonal directions. 
Upon minimizing this action we should obtain the density profile of the vortex dressed 
with the quantum fluctuation of the gauge field. In order to actually minimize this 
action, we must face the problem of solving a system of a non-local and non-linear 
integro-differential equations. This work is in progress and will be presented elsewhere. 

The action Sy should not be confused with the dual actions developed in references 



fjfj and IpTfl . Those dual actions represent vortex densities, whereas the action in the 



present work represents the particle density in the presence of a vortex. In order to 
obtain an action for vortex excitation with our formalism, we would have to introduce 



singularities into configuration space in a way similar to that of ref. p8| . This study 
is under development and it is certainly beyond the scope of this paper. 
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A.l Expansion of the Kinetic term 

In order to deduce the dynamics of density fluctuations, we must expand the kinetic 
term of (|34]) by performing a functional Taylor expansion around the function p(x) = p. 
Here, we would like to sketch the principal steps in developing such expansion. As 
indicated in (|65|) and (|66"D, the kinetic term of (|34]) can be written as 



where 



d 6 z x d 6 z 2 5p(z 1 )K.{z 1 ,z 2 )5p(z 2 ) (92) 



V Z1 • (p{ Zl )V zl K,{zi, Vi)) = S(zt - z 2 ) (93) 



The kernel /C is a functional of p(x) given impicitly by eq. fl9"5|). Putting p(x) = p+Sp(x) 
and expanding up to the first order in Sp, we find 



8K(z x ,z 2 ) 



Sp(x) + ... (94) 



K(z 1 ,z 2 )[5p] = K.(z l7 z 2 )\ + / d 3 x 

J U P^) p=p 

Calculation of the first term is straightforward. From (|93|), we have 

V 21 • (pV Zl K{ Zl ,z 2 )\ p= ^ = pVl K{ Zl ,z 2 )\ p=p = 8( Zl - z 2 ) (95) 

Thus, 

K{z x ,z 2 )\ p=p = \g(z x ,z 2 ) (96) 

where G(z±, z 2 ) is the Green's function of the Laplacian operator. 

The second term of ( |94]) is more involved. Functionally differentiating equation ( |93| ) 
with respect to p we find 

V Z1 ■ (p{zx)V Zl K(z Xl z 2 ) 



5p{x 



= V 21 • (5( Zl -x)V Zl )C(zi,Z2))+V Zl ■ (pizx^^-^Kiz^z^ = (97) 

Observing that the operator V z ■ (p(z)V z ..)j = J dy K,~ l (z,y)... we can rewrite 
equation (^) as 

J dy K,-\ Zl ,y)^-^lC{y,z 2 ) = -V Z1 ■ (5{z x - x)V zx K(z u z 2 )) (98) 
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Inverting this equation we find 
5 



Sp(x) 



K,(z 1 ,z 2 ) = - J dy JC{zi,y) {V y 5(y - x) ■VyJC(y,z 2 ) + 5(y - x)V 2 y K,{y, z 2 )} 



and integrating the 5 — functions we finally obtain 

5 



6p(x 



-K.(z 1} z%) = V x K(zx,x) ■ V x K.(x, z 2 ) 



(99) 



(100) 



This is a functional differential equation for /C. Finding its solutions for arbitrary p(x) 
is a very difficult task, but for our purposes only the special case p = p needs to be 
considered. Using (^) and ( |100p we have 



SK.(z 1 ,z 2 ) 



8p(x) 



-VG{ Zl ,x)-VG{x,z 2 ) 



(101) 



that, of course, coincides with equation (g9| 
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